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ABSTRACT 

An optical model description of momentum transfer in relativistic heavy ion 
collisions, based upon composite particle multiple scattering theory, is presented. 

The imaginary component of the complex momentum transfer, which comes from 
the absorptive part of the optical potential, is identified as the longitudinal 
momentum downshift of the projectile. Predictions of fragment momentum 
distribution observables are made and compared with experimental data. Use of the 
model as a tool for estimating collision impact parameters is discussed. 

1. INTRODUCTION 

Since the pioneering experiments on relativistic heavy ion fragmentation 
using carbon and oxygen beams , 1 ' 2 attention has been directed toward 
understanding the underlying mechanisms of fragmentation processes. Over the 
past two decades, a substantial body of literature has resulted from studies of these 
phenomena, and several excellent reviews have been written . 3 " 6 Perhaps the most 
significant findings of the early experiments were the observations that the 
fragment momentum distributions were Gaussian in the projectile rest frame, and 
that the isotopic production cross sections factored into a product of target and beam- 
fragment terms. Initial attempts to explain these phenomena utilized a statistical 
model to describe the reactions . 7 ' 9 This later evolved into a two-step model called 
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abrasion-ablation 10 where the abrasion stage can be formulated using geometric 10 ' 11 
or quantum mechanical arguments. 12 ' 13 In the present work, we use the impulsive 
excitation energy ideas of Fricke, 14 within the context of composite particle multiple 
scattering theory, to derive a method for predicting momentum transfers occurring 
in relativistic heavy ion collisions. This momentum transfer is a function of impact 
parameter. A new feature of this work is that the momentum transfer is a complex 
quantity. The real component is the usual transverse momentum transfer resulting 
from elastic scattering. The imaginary component is explicitly shown to be the 
longitudinal momentum transfer, or downshift, arising from the absorptive part of 
the complex optical potential. Using this formalism, projectile nucleus fragment 
momentum "downshifts" resulting from the dynamics of the nuclear collision can 
be calculated and compared with laboratory beam measurements. In addition, 
modifications to the widths of the momentum distributions can be estimated using 
the formalism. 

The outline of the paper is as follows. In section 2 the dynamical momentum 
transfer expression is derived, and representative calculations of momentum 
transfer as a function of impact parameter are presented. In section 3 the 
connections between collisional momentum transfer and fragment momentum 
downshifts/ widths are made. A method of choosing appropriate impact parameters 
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for each fragmentation channel is then described. Next, calculations of momentum 
downshifts for fragments produced by oxygen nuclei colliding with various targets 
are made and compared with experimental data. 2 We also compute widths of 
momentum distributions for 139 La fragments and compare with recent experimental 
measurements. 16 In section 4 we propose a method for using the momentum 
transfer model to estimate collision impact parameters. Finally, in section 5 we 
conclude by summarizing the current status of model development and discuss 
future directions for research. 

2. METHOD OF CALCULATION 

In reference 17, a coupled-channels Schrodinger equation for composite 
particle scattering, which relates the entrance channel to all of the excited states of 
the target and projectile, was derived by assuming large incident projectile kinetic 
energies and closure of the accessible eigenstates. The equation is written as 

(v 2 + k 2 ) (x) = 2m A p A T (A p + A T ) * l £ V np> n y (x) y nV (x) 

t ' 

nH 

where the subscripts n and p (primed and unprimed) label the projectile and target 
eigenstates; m is the nucleon mass; A p and A T are the mass numbers of the 
projectile and target; £ is the projectile momentum relative to the center of mass; 
and x is the projectile position vector relative to the target. In terms of the nucleon- 
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nucleon scattering t-matrix t a j, and the internal state vectors of the projectile gP ^ ^ 
and target it was also demonstrated that the potential matrix is expressible as 

Vnjx, n y (x) = (gS gj V opt t Pi i T , x 


P T 
gn'gp' 


( 2 ) 


where 


^opt j^P, ^T, x J — S ^otj 

<*j 


( 3 ) 


This same formalism can be used to investigate heavy ion collision momentum 
transfers. Within the context of eikonal scattering theory, the solution to the 
Schrodinger equation 

H1'(x,£p,4 T ) = E4'(x,£p,£ r ) 
is 



1 

t 

r 

N 

( x > ^P, ^t) = ( 2lt )~ 2 eX P 

* y J V opt ( x » £p. £tJ dz 


( 5 ) 


SK «p> gj (St) <= ik 1 


where v is the velocity. The total momentum of the projectile is then given by the 
matrix element involving the sum of the projectile single-nucleon momentum 
operators as 



'F 


a p - 

E v P , 0 


a= 1 



( 6 ) 


where the subscript P on the gradient operator denotes that the gradient is to be 
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taken with respect to the projectile internal coordinates. Equation (6) actually 
denotes a potential matrix p n ^ n y in analogy with (2). Therefore, substituting (5) 
into (6) yields 


IV, n v 


= ( 6n (^p) gji (£t) 



f A P - ^ 


e -iS 

-i 1 v P , 0 

1 ) 

e iS 


(7) 


where 

S * J:J V opt (x, ^ P , ^ x )dz'. 

— oo 


( 8 ) 


Equation (7) can be further expressed as 


Pnji, n'n' ~ Po + ( 8n 


/ Ap 

-S V P.aS 


A A, 


(9) 


where the momentum before the collision is 


Po=( g£gl 


A P — 

-i I Vp.a 

a = 1 




( 10 ) 


The total momentum transfer to the projectile is then given by 


Q ~ Pn|i, n'n' Pq ” ( gn gj 


_S Vp, a S 


a 


P T 

gn'Sj' 


(ID 


For scattering near the forward directions, the couplings between various excited 
states is small and the off-diagonal elements in Eq. (11) can be neglected; hence, the 
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momentum transfer can be approximated by 


Q-Pei-Po = 


gogl 


v P>a s 



( 12 ) 


In terms of projectile and target number densities, and the constituent-averaged two- 
nucleon transition amplitude, 18 T, Eq. (12) becomes 

Q (b) = -ApA T | d 3 £p pp (^ P ) | d 3 py (£ t ) 


v P )2;(b + ?+4p-i T )^ 


(13) 


where the integration limit in the longitudinal direction has been extended to 
infinity. The momentum transfer in (13) is therefore only a function of the impact 
parameter of the collision. The projectile and target number densities (p p and p T ) 
are normalized to unity as 


j p (x) d 3 x = 1. ^ 

The constituent-averaged two-nucleon transition amplitude is obtained from the 

first-order t-matrix used in our previous studies 13 of nucleus-nucleus collisions as 


t (e, x) = - ( e / m f 2 a (e) [a (e) + i ] [2xcB(e)] ** exp [ - x 2 /2B(e)] 


(15) 


where e is the two-nucleon kinetic energy in their center-of-mass frame, o(e) is the 
nucleon-nucleon total cross section, a(e) is the ratio of the real-to-imaginary part of 
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the forward scattering amplitude, and B(e) is the nucleon-nucleon slope parameter. 
Values for these parameters taken from various compilations are listed in 
reference 18. 

The dynamical momentum transfer to the projectile, given by Eq. (13), results 
from interactions with the target. A new feature, unique to this work, is that it is a 
complex quantity. The real part of the momentum transfer, which comes from the 
real part of the complex optical potential, is the contribution arising from elastic 
scattering. It is purely transverse. The imaginary component, which comes from 
the absorptive part of the complex optical potential, is the longitudinal kinetic 
momentum downshift. To demonstrate this last assertion, we rewrite Eq. (13) 
symbolically as 

Q = (Q R + iQi)b (16) 

where i is V- 1 and £ is the unit vector transverse to the beam direction. If z denotes 
the unit vector in the incident beam direction, then from elementary complex 
analysis 19 we know that 

ib = -z (17) 

since i is an operator which rotates a unit vector counterclockwise through rc /2 


radians. Therefore, the momentum transfer is 
Q = Q R b-Q,z 


( 18 ) 
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which we relabel for clarity as 
Q = Qj.b-Q,,z. 

Note the similarity of this argument to that of complex indices of refraction in 
electromagnetic wave propagation. From Eq. (13), the transverse component is 

= ~ ApAi> | d 3 £p pp (£ P ) | d 3 ^7 p T (£j) V P J^Re t (b + z' + Ijp - £ T ) — 


( 19 ) 


( 20 ) 


and the longitudinal component is 


Qj j = - ApA T j d 3 ^ P p P (^p)j d 3 Px (^ T ) V p Im t (b + z' + - £ T ) 


( 21 ) 


Calculated momentum transfers obtained using equations (20) and (21) are 
displayed in Figure 1 for 16 0 at 2.1 AGeV colliding with a beryllium target. These 
calculations utilize the harmonic well nuclear densities from our previous 
work . 13 ' 18 From the figure, two features are readily apparent. First, the longitudinal 
momentum transfer is larger than the transverse indicating the primarily 
absorptive nature of the nuclear collision at this energy. Second, the predicted 
momentum transfers decrease rapidly with increasing impact parameter. This will 
be a subject of further discussion in subsequent sections of this paper. 
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3. RESULTS 

The collisional momentum transfers computed using the model described in 
the previous section can be related to experimentally-measured, heavy-ion fragment 
momentum downshifts/widths through considerations of energy and momentum 
conservation. As has been pointed out elsewhere, 8 ' 20 a momentum transfer in any 
direction Qj modifies the width hj of the momentum distribution in that direction 

by 



( 22 ) 


and the mean by 

p/-P j+ £Q, (23) 

From the latter, the longitudinal momentum downshift is given by 

4p ii = p lr p ii=f < 2ii <24) 

where Q ( ( is the magnitude of the longitudinal momentum transfer [obtained from 
eq. (21)], F is the fragment mass number, and A is the initial mass number of the 
fragmenting nucleus. Recalling that Q ( ( is a function of impact parameter, an 
appropriate method for choosing it for each fragmentation channel is necessary. 
Recently a semiempirical abrasion-ablation fragmentation model, NUCFRAG, was 
proposed. 21 Although it assumes simple uniform density distributions for the 
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colliding ions, and a zero-range (delta function) interaction, it does include frictional- 
spectator-interactions (FSI) and agrees with experimental cross section data to the 
extent that they agree among themselves. Also, and most importantly for this work, 
it is easily modified to yield the impact parameters for each fragmentation channel. 
Hence, the procedure for evaluating equations (22) and (24) is to extract impact 
parameters from NUCFRAG for each nucleon removal corresponding to AA = 1, 2, 

3, These impact parameters are then inserted into Eqs. (20) and (21) to obtain 

the corresponding momentum transfers for use in evaluating Eqs. (22) and (24). 
Because NUCFRAG uses uniform densities, uniform densities are also used in 
evaluating (20) and (21). In addition, the zero-range interaction in NUCFRAG is 
simulated for numerical integration purposes in (20) and (21) through the use of a 
very narrow Gaussian form for the t-matrix given by eq. (15). This narrow Gaussian 
is the same width for all collision pairs and therefore is not an arbitrarily adjusted 
parameter. 

Representative calculations for momentum downshifts as a function of 
fragment mass number are displayed in Figure 2 for 16 0 projectiles at 2.1 AGeV 
colliding with targets of Be, C, Al, Cu, Ag, and Pb. These momentum downshifts are 
target-averaged using simple arithmetic averaging. For comparison, the target- 
averaged experimental data from reference 2 are also also displayed. For display and 
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comparison purposes, the latter are also averaged over all isotopes contributing to 
each fragment mass number using 


5>iK 




expt 


1 °, 


( 25 ) 


where Oj is the experimental production cross section for the ith fragment isotope. 
Comparing the theoretical estimates to the experimental data, reasonable agreement 
is obtained considering the simplified form of the nuclear fragmentation model 
used in the calculations and the overall sensitivity of the calculated momentum 
transfer to the choice of impact parameter. Improved agreement is expected if 
impact parameters from a fragmentation model using realistic nuclear densities and 
interactions were available. This is especially true for collisions involving lighter 
ions, such as carbon, oxygen, and beryllium, which are poorly represented by simple 


uniform nuclear distributions. 

Figure 3 displays transverse momentum widths as a function of fragment 
mass number for 1.2 AGeV 139 La fragmenting in carbon targets. The experimental 
data are taken from reference 16. Again, impact parameters from NUCFRAG are 
used as inputs into the momentum transfer expressions [Eqs. (20) and (21)]. For 
consistency with the use of these impact parameters, a narrow Gaussian t-matrix 
and uniform nuclear densities were again utilized in the momentum transfer 
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calculations. From Figure 3, it is clear that the agreement is much better than in 
Figure 2 and probably reflects the fact that a uniform nuclear density distribution is a 
more reasonable aproximation for a heavy nucleus-like lanthanum than for light 
nuclei such as oxygen. 

4. ESTIMATING COLLISION IMPACT PARAMETERS 

Thus far in this work, we have used collision impact parameters as inputs 
into a momentum transfer computational model which in turn has yielded 
estimates of heavy ion fragment momentum downshifts /widths for comparison 
with experimental data. However, this procedure can be reversed and the model 
used to estimate collision impact parameters from measured momentum 
downshifts for relativistic collisions. Let F be the fragment mass number with 
measured longitudinal momentum downshift AP , ( produced in a relativistic 
collision between a projectile nucleus (mass number A) and some target. Then, 
from eq. (24), the longitudinal momentum transfer to the projectile from the target 
is 

Q| = A.AP||. (26) 

The collision impact parameter can then be estimated from eq. (21) by computing 
Qi | as a function of impact parameter (e.g., in Figure 1) and using Q, , from eq. (26) 
as the entry. To illustrate, consider a collision involving 2.1 AGeV oxygen colliding 
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with a beryllium target. The calculated momentum transfer using realistic nuclear 
densities are displayed in Figure 1. If the measured (hypothetical) momentum 
downshift for the 14 N fragment is 35 ± 7 MeV/c r then eq. (26) yields a longitudinal 
momentum downshift of 40 ± 8 MeV/c. From Figure 1, the corresponding range of 
impact parameters is 6.1 - 6.4 fm. A similar procedure incorporating measured 
momentum distribution widths and Eqs. (22) and (20) or (21) could also be used to 
estimate collision impact parameters. 

5. CONCLUDING REMARKS 

Beginning with composite particle multiple scattering theory, an optical 
model description of collision momentum transfer in relativistic heavy ion 
collisions was derived. General expressions for transverse and longitudinal 
momentum transfer, which utilize a finite-range two-nucleon interaction and 
relativistic nuclear densities, were presented. The theory was used to estimate 
heavy ion fragment momentum downshifts for relativistic oxygen and transverse 
momentum widths for relativistic lanthanum projectiles. The main new feature of 
this work was the identification of the imaginary component of the momentum 
transfer as the longitudinal collision momentum transfer. Finally, the use of the 
model as a mechanism for estimating collision impact parameters was described. 

The present theory is mainly applicable at intermediate or high energies 
because of the use of eikonal wavefunctions and the impulse approximation. At 
lower energies (below several hundred MeV/nucleon), the validity of straight line 
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trajectories and the assumption of a constant projectile velocity is questionable. 
Therefore, to compare theory with experiment at lower energies 22 revisions to the 
model are necessary. In particular, deceleration corrections to the constant velocity 
assumption are being developed. For incident energies greater than 1 AGeV, first- 
order deceleration corrections are small (< 1 percent). As the incident energy 
decreases, however, the first-order corrections increase significantly (over 50 percent 
at 100 A MeV), indicating that higher-order terms must be included. Work on this 
is in progress and will be reported when completed. 
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Berkeley Laboratory, and Frank Cudnotta of the Environmental Measurements 
Laboratory (USDOE) for useful comments and suggestions. This work was 
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Figure 1: Momentum transfer to the 16 0 projectile, as a function of impact 

parameter, for 2.1 AGeV oxygen colliding with a beryllium target. 


Figure 2: Target-averaged longitudinal momentum downshifts as a function of 

projectile fragment mass number for 2.1 AGeV 16 0 colliding with Be, C, 
Al, Cu, Ag, and Pb targets. The experimental data, taken from 
reference 2, are averaged over isotopes for each fragment mass. 


Figure 3: Transverse momentum widths as a function of fragment mass number 

for 1.2 AGeV 139 La colliding with a carbon target. The experimental data 
are taken from reference 16. 
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